Long-range energy-level interaction in small metallic particles 
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We consider the energy level statistics of non-interacting electrons which diffuse 
in a d- dimensional disordered metallic conductor of characteristic Thouless energy 
E c . We assume that the level distribution can be written as the Gibbs distribution 
of a classical one- dimensional gas of fictitious particles with a pairwise additive in- 
teraction potential /(e). We show that the interaction which is consistent with the 
known correlation function of pairs of energy levels is a logarithmic repulsion for level 
separations e < E c , in agreement with Random Matrix Theory. When e > E c , f(e) 
vanishes as a power law in e/E c with exponents —\,—2, and — | for d = 1, 2, and 
3, respectively. While for d = 1,2 the energy-level interaction is always repulsive, in 
three dimensions there is long-range level attraction after the short-range logarithmic 
repulsion. 
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A statistical description of the hamiltonian H of a complex system is provided by random 
matrix theory. A key feature in this theory is the spectral rigidity of the energy levels: their 
distribution P(E\, E2, ■■■En) formally coincides with the Gibbs distribution of the positions 
of a one-dimensional gas of N classical particles with a repulsive logarithmic interaction, 

P({E n }) = Z- 1 exp {-/3H{{E n })] , (1) 

H{{E n }) = - J £\n\E i -E j \+ J £V(E i ) . (2) 

i<j i 

Here Z is a normalization constant, and V(E) is a confining potential. The parameter j3, 
playing the role of an inverse temperature, depends on the symmetry class of the ensemble 
of random hamiltonians [0] . A method which yields such a distribution consists in assigning 
to the hamiltonian H a probability distribution of maximum information entropy given a 
spectral constraint 0. For instance, a constraint on the expectation value of J2iE 2 yields 
the gaussian ensembles, where V(E) cx E 2 . Other ensembles, characterized by different 
V(E), result from from other spectral constraints (e.g. the averaged level density). All 
these classical ensembles of random matrices have in common an absence of eigenvalue- 
eigenvector correlations and a logarithmic repulsion between pairs of eigenvalues. 

The use of this theory for the study of electronic properties of small metallic particles was 
introduced by Gorkov and Eliashberg |J . Theoretical support for the logarithmic repulsion 
of energy levels came with the work of Efetov Q . Assuming bulk diffusion of the electrons 
by elastic scatterers and using a supersymmetric formalism, he obtained for the spectrum of 
small metallic particles the same correlation function as in classical sets of random matrices. 
Subsequently, addressing the connection between universal conductance fluctuations and 
the universal properties of random matrices, Al'tshuler and Shklovskii || (A&S) showed 
that for energy separations \E — E'\ greater than the Thouless energy E c , the correlation 
function deviates from classical random matrix theory. The Thouless energy E c = Dh/L 2 is 
inversely proportional to the time t erg it takes an electron to diffuse (with diffusion coefficient 
D) across a particle of size L. The results of the diagrammatic perturbation theory of 
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A&S were recently rederived by Argaman, Imry, and Smilansky, using a more intuitive 
semiclassical method [[7]. 

One would not expect that the logarithmic level repulsion in Eq. (0) holds for levels 
which are separated by more than E c . What is then the long-range energy-level interaction 
in small metallic particles? This is the question addressed in this paper. In a sense, this 
is an inverse problem in statistical mechanics: given the pair correlation function, what is 
the interaction potential? Our analysis applies a recently developed functional-derivative 
technique to compute correlation functions in random-matrix ensembles with an arbitrary 
two-body interaction potential ||. The restriction to two-body (i.e. pairwise additive) in- 
teraction is our single assumption. We find that E c characterizes a cross-over between a 
short-range logarithmic repulsion and a novel long-range part which decays as a power law, 
with a dimensionality-dependent exponent. The interaction remains repulsive for dimen- 
sions 1 and 2, but exhibits a long-range attractive part after a short-range repulsion in 3 
dimensions. 

The starting point of our analysis is the Gibbs distribution (|I]) with an arbitrary two- 
body interaction f(\E — E'\) in the fictitious "hamiltonian" TC, 

n({E n }) = J2 f(\Ei - m + £ V{Et) . (3) 

i<j i 

The mean eigenvalue density (p(E)) (J2 n $(E — E n ) (where (. . .) denotes an average with 
weight P({E n })) is related to V(E) and f(E — E') by an integral equation, valid || to the 
leading order of a 1/N expansion: 

/oo 
dE f (p(E f ))f(\E-E'\) + c. (4) 
-oo 

The constant c is to be determined from the normalization condition J dE(p(E)) = N. Eq. 
(|j) has the intuitive "mean-field" interpretation (originally due to Wigner), that the "charge 
density" (p) adjusts itself to the "external potential" V in such a way that the total force 
on any charge E vanishes. Dyson || has evaluated the first correction to Eq. (^), which is 
smaller by a factor iV _1 In N. 
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The density-density correlation function defined by 

K 2 (E, Ef) = (p(E))(p(E')} - (p(E)p(E')} , (5) 
can be expressed functional derivative |§, 

Eq. @ is an exact consequence of Eqs. ([[]) and (Q). Physically, it means that correlations 
between E and E' are important when a modification of the potential at E' has a substantial 
impact on the mean density at E. Combining Eqs. (^) and (||), one can see that K 2 (E, E') = 
K 2 (\E — E'\) is translationally invariant and independent of the confining potential V(E), 
depending on the two-body interaction f(e) only. This property is at the heart of universality 
in random- matrix theory ||,[|. 

By Fourier transforming the convolution (§), the time-dependent two-level form factor 

K 2 (t) = deK 2 (e) exp (7) 

can be written as 

KM = j ww = ~Wm ■ (8) 

This relationship gives us the prescription for obtaining the eigenvalue interaction f(e) from 
the density-density correlation function K 2 (e). 

For disordered systems in the weak-scattering limit, A&S have shown by perturbation 
theory that the correlation function is given by 

K 2 {e) = j^Re £ (e + ihDq 2 + ^)" 2 , (9) 

for energies e large compared to the level spacing A, and small compared with the energy 
scale h/r e , associated with the elastic scattering time r e . The factor s = 2 accounts for the 
spin degeneracy of each level, 7 is a small energy cutoff (to account for inelastic scattering) 
and the parameter f3 equals 1 (2) in the presence (absence) of time-reversal symmetry (/3 = 4 



for time-reversal symmetry with strong spin-orbit scattering) . The sum is over the eigenval- 
ues of the diffusion equation for the sample, assumed to be a d- dimensional parallelepiped 
with sides (q 2 = n 2 Y.fi=i( n ^/ L^) 2 ). In what follows we put s — 1 and 7 = 0, ignoring 
the spin degeneracy and the small-energy cutoff. The Fourier transform of Eq. (|9]) is 

m) = - A £ exp f-^ 2 ^ E(v^) 2 1 • (io) 

The long- and short-range limits K l 2 (t) and iff (i) of the form factor (|T0| ) can be obtained 
in closed form JF]. The cross-over time scale is the so called ergodic time t erg = L 2 /D, which 
it takes an electron to explore the whole available phase space (for simplicity we work with 
a hypercube, = L for all /i). The ergodic time is related to the Thouless energy by 
E c = h/t crg . For times t ^> t crg the first term of Eq. (PH| ) dominates the sum, while for 
t ^ t erg one can convert the sums over n M into gaussian integrals. The resulting long- and 
short-time limits are J7| 

ui Td\ t \-d/2 

w = -jk^w ■ (12) 

For analytical work it is convenient to have an expression which smoothly interpolates 
between these two asymptotic limits. We will use the interpolation formula 

K 2 (t) ~ K l 2 (t) + K s 2 (t) . (13) 

The approximation (|13D differs from the full expression (0) for t ~ t eTg , but it is accurate for 
t either much smaller or much larger then t erg . This is sufficient for the purpose of obtaining 
the asymptotic behavior of the interaction potential. 

Combining Eqs. @ and (|13"D, the interaction potential fd(z) for d dimensions can be 
written as 



f d/2 

fd(z) = I dt cos (at) — i (14) 

/o V ; t{l + t d / 2 ) v ; 



where a = e/ (47rE c ) is the dimensionless energy variable. The integral (|14]) can be evaluated 
numerically for all a, and analytically in the small- and large-a limits. For |a| <C 1 we can 
approximate cos at ~ 1 and cut the upper integration limit at l/|a|, which readily yields 
the short range universal logarithmic interaction fd{s) — — ln|a|. For |a| — > oo the high- 
frequency oscillations of cos at average the integral to zero, in a way which depends on the 
dimensionality. The easiest case is d = 2, where the integral can be evaluated in a closed 
form, 

— ~ sin | a | si|a| — cos (a) ci(a) , (15) 

which behaves as — In \a\ — C for small |a| (C is Euler's constant), and 1/a 2 is the dominant 
term of an asymptotic expansion of Eq. (|T4] ) for \a\ ^> 1. We therefore recover the short- 
range logarithmic repulsion and find that the interaction remains repulsive in the whole 
energy range. For d = 3, the asymptotic limits of the interaction can be obtained by 
considering the auxiliary function 

f 00 sin at ftfe 

h{a) = L dt — TTW' (16) 

which satisfies h'(a) = /3(a) as well as the differential equation 

»// % 1 1 / x 3 f°° 1 sin (aw 2 ) . 
h\a) = -h(a - - / du ^ 3 ; 2 . 17 
a a Jo (1 + u A Y 

The second term on the r.h.s. becomes —1 in the small-a limit, and — 3y/n\2a\^ 3 ^ 2 in the 
large-a limit. We thus obtain /3(e) ~ — In |ot| — C for \a\ 1 and —y/7i\2a\~ 3 ^ 2 for |a| ^> 1. 
Therefore, in d = 3, we have an attractive eigenvalue interaction for large separarations. 
Using a similar procedure for d = 1 we obtain the same short-range logarithmic repulsion 
— In I a I — C, which crosses over to an algebraic repulsion fi(e) ~ y^^al -1 / 2 for a ^> 1. 

In Fig. 1 we compare a numerical integration of Eq. flU]) with the asymptotic expressions 
derived above, which we summarize: 

Z^) = -ln|^|-C ifN«£ c , (18) 



/ \ ,. . _ -1/2 



(f) ,/2 


£ 




4ttE c 



fs{e) = ~\ (f 



► if |e| > E c . (19) 



1/2 



-3/2 



A-kE c 

Fig. 1 shows the crossover from the universal logarithmic short-range repulsion into the 
novel long-range power- law regime (repulsive for d — 1, 2 and attractive for d — 3). We can 
see that for a given E c the departure from the universal regime occurs earlier as we go to 
lower dimensions. 

In summary, we have calculated in the metallic regime the dimensionality-dependent 
long-range part of the energy-level interaction. From this interaction, one could in principle 
calculate n-point correlation functions for arbitrary n. Our method is based on a general 
relation || between the density-density correlation function and this interaction. We use it 
in a particular case where this correlation function is known from diagrammatic perturbation 
theory |J (or an equivalent semiclassical theory Q). The validity of our results is restricted 
to the validity of these perturbative or semiclassical approaches: £ > A, and e <C %/r e . 
Since Efetov has shown that random matrix theory remains valid for e <C A, the universal 
logarithmic repulsion which we recover must be also valid for these small energy separations, 
though either perturbation theory (Eq. ^|), or our method based on an asymptotic large-N 
approximation miss fine structure on the scale of A. 

The condition e <C fi/r e limits the non- universal algebraic decay which we find for 
e > E c = h/t CTg . These non-universal interactions result from the non-ergodic electron 
dynamics for t < t eTg . The non-ergodic dynamics in our problem is unbounded diffusion in 
d dimensions. For times smaller than r e the electron motion is ballistic, a behaviour which 
is not considered in our theory. 

Chaotic billiards with ballistic motion between the boundaries constitute another exam- 
ple for the applicability of random-matrix theory. In the semiclassical limit the energy level 
correlations are also correctly described by the classical ensembles for energy differences 
smaller than the inverse period of the shortest periodic orbits of the system flO| , [TT|j . The 
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essential difference with the small metallic particles considered in this work is that these 
have a well-defined statistical regime of diffusive motion between r e and t erg . 

It is interesting to consider the scale dependence of the interaction potential. In d = 3, 
the level spacing scales as A oc L~ 3 while E c oc L~ 2 , so that A <C E c as L —>■ oo. Therefore, 
if we measure \E — E'\ in units of A, the interaction f(E — E') scales with L towards the 
universal random matrix repulsion for about L nearest neighbour levels for three dimensional 
conductors. However, the total number of levels being proportional to L 3 , the relation (|J) 
between the average density (p(E)) and the confining potential V(E) still differs from the 
usual expression (i.e. with a logarithmic interaction) in the thermodynamic limit. 

Our analysis is restricted to metallic particles which are small compared with the lo- 
calization length. In the thermodynamic limit, electrons are always localized for d = 1 or 
2 (except for (3 = 4 in d = 2 at low disorder). Anderson localization occurs also in three 
dimensions for large disorder. In these cases, our perturbative starting point Eq. @ is 
no longer valid. In the presence of eigenvector localization, f(e) probably scales with the 
system size towards a delta function (uncorrelated levels in the limit of strong localization). 
An interesting issue that we postpone for future studies is to see if the mobility edge is 
characterized by some scale invariant interaction f c (e). 

In this paper we have only considered the eigenvalue statistics. The classical random 
matrix ensembles are invariant under canonical (orthogonal, unitary or symplectic) trans- 
formations. The logarithmic level repulsion is related to the maximum randomness of the 
eigenvectors. In the localized regime, it is clear that the absence of a logarithmic level re- 
pulsion has its physical origin in the localization of the eigenvectors in different parts of the 
system. An important issue would consists in identifying what kind of modification of the 
eigenvector statistics from Porter-Thomas distribution jl| is behind the novel non-universal 
part of the level interaction found in this work. 

Finally, it is likely that studies similar to the one presented here for the hamiltonian 
ensemble will be useful for the ensemble of scattering or transfer matrices [I2~|, and will 
improve our understanding of quantum transport in disordered conductors. 
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FIGURES 

FIG. 1. Interaction potential according to Eq. for various spatial dimensions d (solid), 
together with the short-range logarithmic (dash) and long-range power law (dot-dash) asymptotic 
forms described by Eqs. ( |l8|) and (|l9|). Inset: blow up of the departure from the short-range 
logarithmic interaction. 
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